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Abstract 

In this paper we proposed the alternative test to the two independent and 
normally distributed samples t test based on the cross variance concept. We 
present the simulation results of the power and the error rate of the special 
case of the cross variance which is when the variances of the two samples are 
homogeneous and the t tests. 

The simulation results show that the special case of the cross variance test 
has the power and the error type I rate equal to the t test. This result suggests 
that the proposed test could be used as an alternative to detect whether there is 
difference between means of the two independent normally distributed samples. 
We give some example comparative case studies of the special case of the cross 
variance and the t tests. 
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1 Introduction 


When we have two independent samples and would like to detect whether there is 
difference between means of the two populations then we apply the parametric and 
the non parametric tests. If the samples come from normal distribution, then the t 
test will be the right choice. 

In this paper we will introduce a cross variance concept and present a new approach 
to detect whether there are differences between the means of the two samples based 
on the cross variance. 

This paper is organized as follows. In Section we will describe the cross-variance 
concept and the proposed test, inlcuding the new probability density function. Sec¬ 
tion 1^ will present the simulation results of the power and the error type I of the 
proposed test in comparison to the t test. Some examples will be provided together 
with the t test. Finally, section]^ is the summary and remarks. 


2 An introduction to the cross variance concept and 
the proposed test 


2.1 The cross variance concept 


Definition 2.1. Suppose we have two independent samples, X* and 
and j = 1,2, Their sample mean and variance are denoted by 

Let 




V* = 


2=1 


m — 1 


and V* = 


2=1 


n — 1 


Yj; i = 1,2, ...,m 
X,Y and 14 , 14 . 
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be the cross-variance for each sample X and Y respectively. The cross-variance sam¬ 
ple of groups X and Y is defined as 


T = 



( 1 ) 


where Vf = ^, Vy 



Clearly 


and 


... Er=i(^.-y)^ 

m — 1 


= K + 


m{Y-Xy 


m 


( 2 ) 






2 



n — 1 

n(Y-XY 
n — 1 


Thus T = can be re-written as 


T = - 
2 




-h 

/T.- T.-\0 I 


K, 


14 + K + 


^ ' m—l 

14 


+ 


n{Y-XP 
y ' n-1 

v„ 


214 + 2 = 4 ^ 214 + 2424=42 

^ m—l y n—1 

= Zi -\- Z 2 


(3) 


In what follows, we assume that 

1. the sample sizes are equal 

2. Xi and Yi are i.i.d. normally distributed with unknown means and known 
variances 
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It follows that 


{n-m (n-l)l 4 ,2 „2 


erf 


X(n-l) 


? 9 


X(n-i)> and 




X(i) 


Therefore Equation ([^ can be written as follows 


T = Zi 


(n-l)q, 

^2 


2(!fi:ih4 2 




O'! 


(y-x)^ 


+ 


{n-l)Vy 


> (’^-i)y 

^ 9 


+ 2 


(o'S+o'?) 


( 4 ) 


(o'i+o'2) 


(y-x)^ 


where 


(n-l)Vf, 


^1 = 


2(n-l)14 , 2(, 
0-2 “T 




’■i+o'y) 


(y - x)^ 


u 

2U + 2abV' 


and 


^2 = 


{n-l)Vy 

-I 


s 


2(n—1)V^ I 2{a'^+(T^) ^ 

^ > 


2 2S + 26cy 


with 


and 


_ (n - 1)K ^ _ (n - l)y^ 


erf 


erf 


y = 


n(y -X)^ 


a = —, 6=(a^ + ay), c= — 


erf 


erf 


Hence Equation Q can be written 


as 


T — Z 1 + Z 2 — 


U 


S 


2U + 2abV 2S + 2bcV 


( 5 ) 


To compute the distribution of T in Equation ([^, it can be done by considering the 
fact that 


1. U,V and S are independent 
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2. Zx and Z 2 are dependent 


In this paper we will describe the first case, where we consider that f7, V and S are 
independent. 


2.2 The proposed test 


Under normality asumption of X and V then [/, V and S are independent, where V 
is X(i) distribnted and U,S are X(n-i) distributed. From the equation (jS), suppose 
V = Z3 from here we get that U = ^ ^ Jacobian of this 

transformation is 


1^1 = 


dU 

dU 

dU 


2abZz 

0 

2a6Zi 


dZi 

dZ2 

dZs 


(1-2Zi)2 

(l-2Zi) 


dS 

dS 

dS 

_ 

0 

2bcZ'^ 

2bcZ2 


dZi 

dZ2 

dZs 


(1-2Z2)" 

(1-2Z2) 

(d 

dV 

dZi 

dV 

dZ2 

dV 

dZa 


0 

0 

1 


Aab^cZ^ 


((1-2Zi)(1-2Z2))' 


The joint probability density function of Zi, Z 2 , Z^ is 


/zi,Z 2 ,^ 3 (U, ^ 2 , ^3) 

[ 2abziZz\ [ 2bcz2Z3\ 


(4a6^c) 


n — 1 
2 






2n-l 
/ 2 


-1 

^3^2^ l-2zi^ 1-2 z2 J 


■ir(i)r(V) 


n+1 n+1 

(l-2zi)—(1-2 z2) — 


( 6 ) 


The joint density function of Zi, Z 2 is the marginal probability function of Zi,Z 2 
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from the equation ([^ above. It is computed as follows: 


/zi,Z2(^U ^ 2 ) 


CXD 

j fZ t_,Z2,Z3{zi, Z 2 , Z3)dz^ 

0 


(4a6^c) 


n —1 
2 


n —1 


-1 


n — 1 


-1 


2n-l 
/ 2 


n+1 


2”Ur (i) r (ti)" (1 - 2.-0“ (1 - 


-1 

-23^2^1-221 


bc 22 N 

^-'^^ 2 ) dz^ 


n — 1 n — L 1 

(4a6^c) ^ T {n — Zi^ z. 




-u ( 1 ) r (1 - 2z,)’^ (1 - 2z,)’^ (i + (i?!So + (Ss)’ 


( 7 ) 


Therefore the cumulative distribution function (cdf) of T <ts computed as follows 


Frit) 

= P{T < t) 


00 t—Zi 



fz3,Z2{Zi,Z2)dZ2dZi 


—00 —00 

t t — Zi 



n — 1 _^ n—1 _^ 


(4a6^c) ^ r (n — I) 2]^ ^ 22 ^ dz2dzi 


f f ^ O n-j-l / 

i i 2 ”-ir(yr(!!fi) ((i- 2 z,)(i- 2 z,)) — [\ + 


n— 

ahz\ I hcz2 * 


( 1 - 22 : 1 ) ^ ( 1 - 222 ) 


t-i 


n —1 TV 

(4a6^c) ^ r (n — I) 2 ]^ 

2-ir(i)r(^)'( 1 - 221 )^ 


'Qdzi 


( 8 ) 


where 


B 



n — 1 

^ 2 ^ 


-1 


(1 


n+1 

22o') 2 (i2o 



ab2:i I 

( 1 - 221 ) 


bcz2 

( 1 - 222 ) 



( 9 ) 


To compute the integral at Equation ([^, hrst we will simplify this 

/1 I abzi I bcz2 \ ^ ^ pQ /^1 2(fec-(l-2(l-afe-bc)2i)) ^ \ f l-2{l-ab)zi 

\2 (1-221) (1-222) J l-2(l-a6)2i \ 2(1-221) 


. Therefore the 
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Equation (|^ becomes 


(2(1-2z,)) 


”-3 




n+1 


(1 - 2 Z 2 ) 2 dZ2 


(1 — 2(1 — ab)ziy 


2 J /. 2 {bc-(l- 2 {l-ab-bc)zi)) ^ 

0 1 1 + l- 2 {l-ab)zx ^2 


( 10 ) 


t — Zi 


-1. . ^+1 


The integral J 


22 ^ (1-222) 2 Ci 22 


/T I 2(!)c-(l-2(l-ai)-be)zx)) _ 
l-2(l-a6)zi 


— is written as 


t — Z\ 


n —1 

^2 


»-■(!- 2.-3)'^ (1 + 2(6c-(l-^2(l-a6-6c).-.)) ^^ 


1 — 2(1 — a6)^i 




dz2 (11) 


By considering the binomial expansion then Equation ([TT| can be represented 


as 


00 00 

k=0 1=0 

t — Zi 


^^^^{bc-l + 2(1 -ab- 6c)^i)' 


(1 — 2(1 — ab)zi) 


X 


-+k+l 


dzo 


00 00 


EE 

k=0 1=0 


2 k+i^_^^k+i ^^^^{bc-l + 2(1 -ab- 6c)xi)' 


(1 — 2(1 — ab)zi) 


X 


^^+k+l 


( 12 ) 


6 















Therefore, 


B = 


(2(1-2ji))’ 


(l — 2(1 — ab)ziy 


EE 

. k=0 1=0 _ 

^ + k + l 


k 


-X 


n+l 


n 


i + / 




i j 


{be — 1 + 2(1 — ab — bc)ziy 


(l — 2(1 — 06)2:1)* 


■X 


(13) 


and 


Frit) = 


n — 1 

(Aab'^c) 2 r (n - I) 

r(yr(iifi)^ 

(-2)'I I )l 


X 


EE 

k=0 1=0 


m-r 


{—2{bc — l)yt "2 


^^ + k + l 


G 


(14) 


where 


G = 


n — 1 

2^1 




n+l 




(l-2zi)t-i (l-2(l-a6)xi) 




X 


2(1 — a6 — bc)zi , , 

bc-l I 


(15) 


Again, by considering the binomial expansion then Equation (15) can be written as 
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follows 


OO OO OO 


m=0 p=0 q=0 L 





^2+m+p+q ^ _ ^}_ 


OO OO OO 







m=0 p=0 q=0 L 

B {^ + m + p + + k + l) 

^^^^+m+p+(j 


2),„/ 2(l-«6-6c) y(2(i^^)),^ 


6c — 1 


( 16 ) 


n — X 

^ ^ {^ab^c)~T{n-l) 

r(h) 



^^y^^+k+i 


^ + k + l 



X 



^-^ + k + l 

(_2)m 

\ be — ] 


-2{bc - 1))V+^- 


X 


■m—p—q 


X 


Furthermore, from the Equation ([T^ then the probability density function (pdf) of 
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T is 


... k )[ I , ^ 

' r /'i'\ r 2-^ 2-^ Z-^ Z-^ 2-^ " 

■^V2/'^V 2 / fc=0 1=0 m=0 p=0 g=0 


^ + k + l 


{k + l-m-p- q)t’^+^-^-P-‘i{-2{bc - 1))^ 0 C ^ ^ g ^ ^) 

2(1 — ab — bc)\^, , ,,,„ fn — 1 n + 1 , \ , , 

-^-j {-2{l-ab)yB (^— + m + p + q,^— + k + lj (18) 


The pdf of T also can be computed as follows: 

t 

frit) = j fz^,Z2{.zi,t - zi)dzi 

— OO 

t 

= j fZuZ^it - Z2, Z2)dZ2 

— OO 


(19) 


Now we already have the pdf and cdf of T from where we can compute the T statistics 
value for the hypothesis testing, that is the hypothesis null about the equality of mean 
of two-groups independent samples is rejected if T < to or P (^T < to) = Frido) < «• 


The computation of Frito) by using the Equation (18) is involving the hve summation 
and hence it is not quite simple. The computation gets simple in the case = a^, 
which will be described next. 


2.3 Special case of the proposed test 

In the case of we estimate the 14 and I4 by the least square estimator of 

the pooled variance = hk±4_ 

If we use the least square estimate as the estimator of 14 and I4 therefore the Equa- 
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tions ([^ and (|^ becomes 


T* = - ^_, 

V^^+Vy , 

2 ^ n-1 

u* 

~ U* + AV* 


where U* 



and 1/* 


n(r-x)^ 

20-2 


(20a) 

(20b) 


The pdf of T* is derived from the ratio of linear combination of chi-square random 
variables [1]. First, let y = 1 + 4^, where V* is distributed xfi) and U* is distributed 
X 2 {n-i)- Second, the pdf of T* is computed by taking T* = y. 


In the folowing computation, the chi-square distribution will be represented as the 
Gamma distribution. Therefore, we have that V* is Gamma distributed with param¬ 
eters ^ and Pi = 2. U* is Gamma distributed with parameters a 2 = {n — 1) 

and /32 = 2. U* and V* are independents. 


Suppose G = ^ and if we take U* = H, then we get V* = GH. Further we got 
the Jacobian of this transformation variable random is h. Because V* and U* are 
independents, then the joint probability function of G and FT is 


fG,H{g,h) — fv*,u*{gh,h).h 


( 21 ) 


where 


fv*,u*{gh,h) — fv*{gh).fu*{h), 


fv* {gh) 


fu*{h) 




gh 


PTV{ai) 

(h)“2-ie"^ 


and 
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Therefore 


fG,H{9, h) 


/3“T(ai) ■ /3rr(a2) 
/3«i+“2r(ai)r(a2) ’ 


h 


/Si 


1^2 — /3 


We want to determine the pdf of 5 / then 


( 22 ) 


00 



fG,H{g,h)dh 


B{ai,a2) (1 +5f)“i+“2 




(23) 


G is distributed beta of second kind. 


The next step is determining the distribution of Y. We dehne y= 1 + AG and by 
using the transformation random variable, where G = the pdf of Y is computed 
as follow 




(24) 


The pdf of T* = y is obtained 




B{ln-1) {3+Y)n-ht*^ 

J- - -r, o<r<i 

B [^,n- 1 ) (1 + 3f*)” 2 


(25) 
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Furthermore the cdf of T* analtically is computed as 


Ft*{Q = j fT*{f)df 

0 

t* 

gn-l 

“ J (l + 3t*r-i 


4n-l 

B (1,71-1) 

4n-l 





^ — 1 + k, 



X 


pbeta 


^to,n - 1 + A;, 



(26) 


We reject the null hypothesis of the equality of mean of two independent samples if 
T* < oi Pit* < T*) = FT*iT*) = p-value < P (t* < = a. 

Observing that (2(n — 1)) is the square of a random variable having f 2 (n-i) 

distribution, a simple calculation shows that the same holds for the random variable 


J = 



(27) 


This statistic J can also be used to test the hypothesis = /iy and the critical values 
can be computed from the t table. It also follows that the J has a limiting normal 
distribution as n —)■ cx) 


3 Simulation study 

In this section we will describe the results from the two simulation studies of the t and 
the proposed tests. First simulation is done to measure the power of the proposed test 
is conducted at iV = 500 times and a = 0.01. The results are presented at Section 
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3.1 In this simulation, we consider various possibilities regarding the sample sizes 


and variances. The results of the simulation are divided into groups according to 


1 . variance: low (S' = 0.2), medium (S' = 1.2) and high (S' = 7), 

2. sample size: low (5 and 25), medium (100) and high (500) 


Section |3.2| describes the second simulation results on measuring the rejection’s rate 
under the null hypothesis of the t and the proposed tests, based on the N = 500 
times simulation and a = 0.01. In the simulation, we choose = hv = h = 9-2 and 
ax = cty = cr are three values of variance 1.25, 3.5 and 10 which represents the low, 
medium and high variance respectively. 


Section 3.3 gives some examples on how the proposed test works on the data. 


3.1 Power of the test 

The simulation study to compute the power of the proposed and the t tests is con¬ 
ducted as follows: 


1. Choose the = cry = cr of the two-groups independent samples 

2. The simulation under null hypothesis 


(a) Generate n random samples normally distributed with mean and standard 
deviation, Hxq and cr 

(b) Generate n random samples normally distributed with mean and standard 
deviation, hyq and cr 


(c) Gompute their mean and variance samples 


(d) Gompute values, based on the equation (20a). 


3. The simulation under alternative hypothesis 
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(a) Generate n random samples normally distribnted with mean and standard 
deviation, and a 

(b) Generate n random samples normally distribnted with mean and standard 
deviation, /iy^ and a 


(c) Gompute their mean and variance samples 


(d) Gompute t\ values, based on the equation (20a). 


4. 

5. 


Repeat steps (2) - (3), M times 


Gompute the a quantile of t^. Note that fg a can be computed by using 
the quantile of pdf of T* in the Equation (21). 


6 . Gompute the power test of the proposed method = 1- 

7. Gompute the power of t test from samples Xi and Yi 

8 . Gompare the results from steps (6) and (7) 

9. Do steps (l)-(8), for different value of mean and variance 


Figures 3.1, 3.2, 3.3 and 3.4 show the power of the special case of the cross variance 
and t tests, based on the simulation study. They show that the proposed and the 
t tests have an equal power. In the computation, the power of the proposed test is 
provided by implementing the the empirical approach, therefore the results is not as 
smooth as the t test. 
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Figure 3.1: Graphical power both the t and proposed tests, n=5 



Figure 3.2: Graphical power both the t and proposed tests, n=25 


a^4.01,ip100,»i 






Figure 3.3: Graphical power both the t and proposed tests, n=100 


3.2 Error type I rate 


The simulation study to compute the error type 1 rate is conducted as follows: 
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Figure 3.4: Graphical power both the t and proposed tests, n=500 

1 . Choose the n, Mi,a, = IJ^y = and ax = cry = cr of the two-groups 
independent samples 

2. Repeat Mi times 

(a) Compute p-value of test 

(b) Compute p-value of t test 

3. Compute the proportion of t test p-value < a in Mi results 

4. Compute the proportion of the proposed test p-value < a in Mi results 

5. Compare the result of steps 3 and 4 


In Table [T| it is shown the results of the simulation where the error type I rate of the 
proposed and the t tests are equal. 


The distribution of p-values from the proposed and t tests are can be seen in the 


Figures 3.5 3.6 and 3.7 


The next section will give some examples of data analysis concerning the proposed 
test and compare the results with the t test. We used 14 artificial data sets which 
the first 10 were taken randomly from the internet. The data can be seen in Table 
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Table 1; Error type I rate under 500 simulation for the proposed and t tests 


Sample size Variance 

proposed test 

t test 

0.05 

0.01 

0.05 

0.01 


low 

0.056 

0.012 

0.056 

0.012 

5 

medium 

0.062 

0.016 

0.062 

0.016 


high 

0.056 

0.020 

0.056 

0.020 


low 

0.046 

0.010 

0.046 

0.010 

25 

medium 

0.044 

0.012 

0.044 

0.012 


high 

0.038 

0.010 

0.038 

0.010 


low 

0.058 

0.006 

0.058 

0.006 

100 

medium 

0.050 

0.010 

0.050 

0.010 


high 

0.062 

0.012 

0.062 

0.012 


low 

0.038 

0.004 

0.038 

0.004 

500 

medium 

0.038 

0.002 

0.038 

0.002 


high 

0.050 

0.010 

0.050 

0.010 


cros$ variance test ttest cross variance test ttest 



(a) n=5 


(b) n=500 


Figure 3.5: P-values distribution of the proposed and t tests, small variance 



(a) n=5 (b) n=500 

Figure 3.6: P-values distribution of the proposed and t tests, medium variance 

3.3 Some examples 

In this section, some examples data sets are provided and used as an example how to 
take the decision by using the special case of the cross variance test T* (= the T test 
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(a) n=5 


(b) n=500 


Figure 3.7: P-values distribution of the proposed and t tests, high variance 

when ax = o'y)- For this T* test, the first step is making sure the variances of the 
samples are equal. There are some tests for this purpose, in this paper we implement 
the F.test in R. 

The data sets are shown as in Table and the computation results for the F.test are 
shown in Table Tables]^ and [^provide the results of the decision for both methods: 
the proposed and the t tests. 

In the proposed test, the null hypothesis is rejected if t*o < ta or if the p-value of the 
special case of the cross variance test less than a, where t* is the statistics from the 
observed sample and p-value = P (tl < ). In the computation we use a = 0.01. 

From Table 1^ we can infer that data set 13 has unequal variances, therefore it will be 
excluded from further computation. 

From the Table we can see that the p-values and decisions from the proposed and 
the t tests are equal, except for the data set 4. It is the only one data set which has 
different sample size. When the sample size of two samples is different, we provide 
two options of n: the max(ni,n 2 ) or the average between (ni,n 2 ). The example of 
this is shown in Table [H 
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4 Summary and Remarks 


In this paper we have introduced the cross-variance concept, a new test based on 
the cross variance, the special case of the cross variance (as the proposed test) and 
the new probability density functions. A simulation study has been conducted to 
compute the power and the error type 1 rate of the proposed and the t tests. 

The simulation study of the power shows that the proposed and the t tests has the 
same power. Furthermore the p-value and the error type I rate under the null hy¬ 
pothesis of the proposed and t tests are exactly equal. These results suggest that the 
proposed test could be used as an alternative test to detect whether there are differ¬ 
ence between the means of two independent normal populations, where the variance 
and sample size is equal. In case the sample size is unequal we proposed to choose 
n = max (ni, n 2 ) or the average of (ni, n 2 ). 
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Table 2; Data sets, their mean and variance 


Data 

Data 

Mean 

Variance 

1 

X=(5,7,5,3,5,3,3,9) 

Y=(8,l,4,6,6,4,l,2) 

5.000 

4.000 

4.571 

6.571 

2 

X=(0.72,0.68,0.69,0.66,0.57,0.66,0.70,0.63,0.71,0.73) 
¥=(0.71,0.83,0.89,0.57,0.68,0.74,0.75,0.67,0.80,0.78) 

0.675 

0.742 

0.002 

0.008 

3 

X=(42,45,40,37,41,41,48,50,45,46) 

¥=(43,51,56,40,32,54,51,55,50,48) 

43.500 

48.000 

15.833 

57.330 

4 

X=(33,31,34,38,32,28) 

¥=(35,42,43,41) 

32.667 

40.250 

11.067 

12.917 

5 

X=(35,40,12,15,21,14,46,10,28,48,16,30,32,48,31,22, 12,39,19,25) 
¥=(2,27,38,31,1,19,1,34,3,1,2,3,2,1,2,1,3,29,37,2) 

27.150 

11.950 

156.450 

213.524 

6 

X=(26 21,22,26,19,22,26,25,24,21,23,23,18,29,22) 

¥=(18, 23, 21,20,20,29,20,16,20,26,21,25,17,18,19) 

23.133 

20.867 

8.552 

12.552 

7 

X=(520,460,500,470) 

Y=(230,270,250,280) 

487.500 

257.500 

758.333 

491.667 

8 

X=(3,0,6,7,4,3,2,l,4) 

Y=(5,l,5,7,10,9,7,ll,8) 

3.333 

7.000 

5.000 

9.250 

9 

X=(16,20,21,22,23,22,27,25,27,28) 

¥=(19,22,24,24,25,25,26,26,28,32) 

23.100 

25.100 

13.878 

11.878 

10 

X=(91,87,99,77,88,91) 

Y=(101,110,103,93,99,104) 

88.833 

101.667 

51.367 

31.867 

11 

X=(10.11,7.36,6.34,11.83,8.61) 

¥=(3.28, 6.52,2.28,6.66,4.55) 

8.850 

4.658 

4.761 

3.760 

12 

X=(4.79,4.95,2.52,4.98,4.99) 

¥=(7.90,7.51,6.62,7.57,7.49) 

4.446 

7.418 

1.166 

0.227 

13 

X=(3.99,3.98,4.03,4.06, 3.84) 

¥=(6.68,6.25,6.97,5.75,4.01)) 

3.980 

5.932 

0.007 

1.366 

14 

X=(10.16,8.26,16.23,1.44,0.66) 

Y=c(28.06,8.52,25.39,15.45,16.03) 

7.350 

18.690 

41.816 

63.422 


Table 3; F.test decision 


Samples 

Decision 

Samples 

Decision 

Data 1 

equal variance 

Data 8 

equal variance 

Data 2 

equal variance 

Data 9 

equal variance 

Data 3 

equal variance 

Data 10 

equal variance 

Data 4 

equal variance 

Data 11 

equal variance 

Data 5 

equal variance 

Data 12 

equal variance 

Data 6 

equal variance 

Data 13 

not equal variance 

Data 7 

equal variance 

Data 14 

equal variance 
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Table 4: P-values and decision from the proposed and t tests 


Sample 

t test 

proposed test 

p-value 

Decision 

p-value 

Decision 

Data 1 

0.411 

Accept 

0.411 

Accept 

Data 2 

0.054 

Accept 

0.054 

Accept 

Data 3 

0.114 

Accept 

0.113 

Accept 

Data 4 

0.009 

Reject 



Data 5 

0.001 

Reject 

0.001 

Reject 

Data 6 

0.067 

Accept 

0.066 

Accept 

Data 7 

0.000 

Reject 

0.000 

Reject 

Data 8 

0.010 

Accept 

0.010 

Accept 

Data 9 

0.229 

Accept 

0.229 

Accept 

Data 10 

0.006 

Reject 

0.006 

Reject 

Data 11 

0.012 

Accept 

0.012 

Accept 

Data 12 

0.000 

Reject 

0.000 

Reject 

Data 14 

0.039 

Accept 

0.039 

Accept 


Table 5; P-values and decision from the proposed test, Data set 4 


Choice of n 

Least Square 

p-value Decision 

Min 

0.021 Accept 

Max 

0.004 Reject 

Average 

0.009 Reject 
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